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Abstract. In the present paper, integrable semi-discrete and fully discrete analogues of a 
coupled short pulse (CSP) equation are constructed. The key of the construction is the bilinear 
forms and determinant structure of solutions of the CSP equation. We also construct N- 
soliton solutions for the semi-discrete and fully discrete analogues of the CSP equations in 
the form of Casorati determinant. In the continuous limit, we show that the fully discrete CSP 
equation converges to the semi-discrete CSP equation, then further to the continuous CSP 
equation. Moreover, the integrable semi-discretization of the CSP equation is used as a self- 
adaptive moving mesh method for numerical simulations. The numerical results agree with 
the analytical results very well. 
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1. Introduction 

The short pulse (SP) equation 

1 , 

Uxt = u + -{u) xx (1.1) 

o 

was derived by Schafer and Wayne to describe the propagation of ultra-short optical pulses 
in nonlinear media |0Q [2). Here, u — u(x.t) represents the magnitude of the electric field, 
while the subscript t and x denote partial differentiations. The SP equation represents an 
alternative approach in contrast with the slowly varying envelope approximation which leads 
to the nonlinear Schodinger (NLS) equation. As the pulse duration shortens, the NLS equation 
becomes less accurate, while the SP equation provides an increasingly better approximation 
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to the corresponding solution of the Maxwell equations With the rapid progress of ultra- 
short optical pulse techniques, it is expected that the SP equation and its multi-component 
generalization will play more and more important roles in applications. 

The SP equation has been shown to be completely integrable0[4l[5l. The loop-soliton 
solutions as well as smooth-soliton solutions of the SP solution were found in 0[7l(8l. Multi - 
soliton solutions including multi-loop and multi-breathers ones were given in flU. Periodic 
solutions to the SP equation were discussed in IfTOll . 

Similar to the case of the NLS equation [fTTl . it is necessary to consider its two- 
component or multi-component generalizations of the SP equation for describing the effects 
of polarization or anisotropy. As a matter of fact, several integrable coupled short pulse have 
been proposed in the literature [121, 13, jT4[ 15, T611T71 . A complex version of the integrable 
coupled short pulse equation in [TT4l IT5 1 is studied in llT8ll . The bi-Hamiltonian structures for 
the above two-component SP equations were obtained in lfl9ll . 

Integrable discretizations of soliton equations have received considerable attention 
recently [I201I2T1 221 [23J. Integrable semi- and fully discretizations of the SP equation were 
constructed via Hirota’s bilinear method Il24ll . The same discretizations were reconstructed 
from the point view of geometry in lf25l . Most recently, an integrable discretization for a 
coupled SP equation proposed in [ 1.4. 15] was constructed lf26l . 

In the present paper, we consider integrable discretizations of another coupled short pulse 
(CSP) equation proposed by one of the authors lfT6l 


1 , 3 , 1 , 

Uxt = u + -(u ) xx + -v~u xx , 
O Z 

1 , 3 , 1 2 

Vxt = V + - {V^xx + —U V xx . 
o Z 


( 1 . 2 ) 

(1.3) 


It was shown in llT6l that Eqs. (11.21) and (11.31) can be derived from bilinear equations 


DsDyf-f= l -(f 2 -,f 2 ), 

(1.4) 

DsD y f-f= l -(f 2 ~f 2 ), 

(1.5) 

D S Dyg'g= ^{g 1 -g 2 ); 

D S Dyg ■ g = ^(g 2 ~ g 2 ) , 

(1.6) 

(1.7) 

through a hodograph transformation 


x = y— (ln(FF))^ , t = s, 

(1.8) 

and dependent variable transformations 


“ =i ('4)y v=i ( ta I)/ 

(1.9) 


where F = fg. G = fg, F and G stand for the complex conjugate of F and G, respectively. 
Meanwhile, A-soliton solutions of CSP equation in a parametric form are given, and the 
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properties of one-soliton, soliton-breather solutions are investigated in detail in Ifl6l . The bi- 
Hamiltonian structure of CSP equation (ll.2l) - (U.3l) is derived by Brunelli and Sakovich [[T9l . 

The rest of the present paper is organized as follows. In section 2, we propose an 
integrable semi-discrete analogue of the CSP equation via constructing a Backlund transform 
to the bilinear equations of the CSP equation. Meanwhile, A-soliton solution is provided 
in terms of Casorati determinant form. In section 3, starting from two sets of Backlund 
transforms to the bilinear equations of the CSP equation, the fully discrete analogue of 
the CSP equation is proposed by introducing two auxiliary variables. Moreover, /V-soliton 
solution is presented to confirm the integrability. Section 4 is contributed to the self-adaptive 
moving mesh method by applying the semi-implicit Euler scheme to the semi-discrete CSP. 
The paper is concluded by section 5. Appendix A, B and C present the proofs of Proposition 
1, Theorem 1 and Theorem 2, respectively. 


2. Integrable semi-discretization 


We start with two sets of bilinear equations for the semi-discrete two-dimensional Toda-lattice 
(2DTL) equations with the same discrete parameter a 


Id 


x-i 


1 ) X„(fc+ 1) • X n (k) +X n+ 1 (k + l)x„_i(k) — 0, 


~ D x-\ - ^n(k+ \)-x' n (k)-\-x' n+l {k-\-\)x' n _ x {k) — 0, 


( 2 . 1 ) 
( 2 . 2 ) 

which is linked by a Backlund transformation [|27l 

(D x _ t - \)x n {k) ■x' n (k)+x n+ i(k)'t' n _ l {k) = 0. (2.3) 

Proposition 1 The bilinear equations (12. 1 1) — (12.3b admit the following determinant solutions 


x n (x~i,k) = 


x' n (x-l,k) = 


>n\k) 

O*) 


)^\k) 

O*) 

■■■ ^hIn-1 

?°(*) 

♦£?,(*) 

::: e.t 

in\k) 

VnliW 

••• vi +N- 

$\k) 

vi+iW 

••• vi+tf- 

f?\k) 

vi+iM 

••• vi+iv- 


(2.4) 


(2.5) 


where 


with 


tyn\k) = p ?(1 - apt ) k e t ” + q'-( 1 - aqf) k e r]i , 

Vn\k) = p n i{\~Pi){\ -api)~ k e^+ q n i (\-q i )(\-aq i )- k e ! ' i , 

*~>i — Pi X— 1 + ^/0 5 fi; — qi -X— 1 T t|;0 • 


( 2 . 6 ) 

(2.7) 
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Here pi, q t , c,,o and r|,o are arbitrary parameters which can take either real or complex values. 

The proof is presented in Appendix A. 

Applying a 2-reduction condition q t — —pu then we could have each of the x sequences 
become a sequence of period 2, i.e., x„ =c= l n + 2 , f, - x' 2 . Here =c= means two x functions 
are equivalent up to a constant multiple. Furthermore, by choosing particular values in phase 
constants, we can make x„ and x n+ i complex conjugate to each other. Based on the bilinear 
equations with 2-reduction, we construct semi-discrete analogue of the CSP equation by the 
following theorem: 

Theorem 1 The following equations constitute an integrable semi-discretization of the 
coupled short pulse equation (1 1 .2b — (II .3b 


Y s Ofc+t ~u k ) = ^8 k(uk + i+u k ) - ^-(Mjfc+i -u k )(vl +l - v*), 


-J- (Vjt+i -v*) = ^8^(v fe+1 + V*) 
ds 2 

1 / 2 2,2 2 ^ 
-JJ - ~2 y U k+l~ U k + V k+l ~ V k) ■ 


1 

25^ 


(v k+1 -v k )(ul +1 -u 2 k ), 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


Furthermore, A-soliton solution to the above semi-discrete CSP equation is of the following 
form 

' fkg k 


Uk — iln 


, Vk — iln 


Jkgk 

x k = 2 ka - (in (fkfkgkgk)) s , 

8fc = %k +1 %k = 


f fkgk \ 
\fkgk )/ 


( 2 . 11 ) 

( 2 . 12 ) 


+ : 


where j\, gk, fk and gk arc tau-functions defined by 




■ + 


fk XoL, k) , fk H ( , kj , gk x c 


2 ’ 


with 




('*+/'-1) 


1 <i,j<N 


k = 




gk+lgk 


(2.13) 

gk+lgk 

J’ 

gk = xi 

(5-*) ■ 

(2.14) 

1 <i 

J<N ’ 

(2.15) 


$\k) = p^l-api)-^^ + (- Pi ) n ( 1 +a Pi ) 
-k n fir s +^io 




-k~frs+r\K) 


= P?(l -pi)(l-ap i )- k t^’ +i “ + (-pi)"(l+p,)(l +ap,)- t e-& + ''». 

The proof is presented in Appendix B. In the process of the proof, multi-soliton solution 
expressed in determinant form is obvious. Next, we show that the semi-discrete CSP equation 
converges to the CSP equation in the continuous limit. 

In the continuous limit a —> 0 (8^ —> 0), we have 

Uk+l + Uk 


W£_|_i — Uk du 

8 k dx ’ 


—y u ? 


2 


(2.16) 
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dx dxo ^ dbj 
ds ds ds 


k-l 




>+i + v ^+i — — v ^) —^ + v )• (2.17) 


7=0 


Thus 

dr 1 

3, = 3, + ^3 V 4i,--(» ! + v 2 )3,. (2.18) 

Consequently, Eq. (IB.25I) converges to 



u x — u{\ +ul ), 


which is nothing but the first equation of coupled short pulse equation (11.21) . 

It can be shown in the same way that Eq. (IB .261) converges to Eq. (11.31) . the second 
equation of coupled short pulse equation. 


3. Fully discretizations of the coupled short pulse equation 


To construct a fully discrete analogue of the CSP equation, we introduce one more discrete 
variable l which corresponds to the discrete time variable. 

It has been known in 1(291 that the x-functions 




1 <i,j<N 


(3-D 


with 


.(Or t n_„nn A s+ ^ +q >;(]- aqi )-k e^l s+r]i0 


% (k,l) = Pi(l-api)' 


Pi 


Wn\k,l) = p?(l-pi)(l-api) k ( 1 -b —^ e 2 Pi s+ ^° 


Pi 


+q l l \l-q i )(l-aq i )- k [l-b-) e^ +m 


Qi 


satisfy bilinear equations 
'2 


—Ds — 1 Jn(k+ 1,/) -T n(k,l ) +X„+i(k+ 1, l)i n -i {k, l) — 0, 

D s — 1^ x' n (k + 1,1) -x' n {k,l) + x' n+1 (k + 1 ,l)x' n _i(k,l) — 0, 


and 


(2 bD s — 1 )x„(k, l + 1) • x„+i (k,l) + x n (k, l)x n + 1 (k,l + 1) — 0. 
(2bD s - l)x’ n (k,l+ 1) -x' n+l (k,l) + <(M)< +1 (M+ 1) - 0. 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


Here n,k,l are integers, a, b are real numbers, p\, qt, q,o-B/o are arbitrary complex numbers. 












Integrable discretizations and self-adaptive moving mesh method for a CSP equation 


6 


By applying a 2-reduction condition: qi — —p,, we have x n o x n+ 2 , x' n =c= x' ;+9 . We can 
further have 


A,/ =T o(M)» fkj = , gkj = x' 0 (k,l), gk,i = x\(k,l), 

by adjusting phases in tyf l \k,l) and \\f\”\k,l). Here / and g represent complex conjugate 
functions of / and g, respectively. A fully discrete CSP equation can be constructed as 
follows: 

Theorem 2 The fully discrete analogue of the CSP equation (ll.2lHll.3D is of the form 

-^( u k+l,l+l~ u k+l,l~ u k,l+l+ u k,l+ v k+l,l+l~ v k+l,l~ v k,l+l+ v k,l) 

= ( yk+ 1,1 —yk,i+ 1) {uk+ 1,/+1 + u k,i+vk+i ,i +1 +vk,i ) 

+ (yk+l,l+l—yk,l)( u k,l+l+ u k+l,l+ v k,l+l+ v k+l,l) , (3-6) 

= (Zfc+1,/—Vjfc+l^+t —Vjfc,/) 

+ (Zk+l,l+l-Zk,l)( u k,l+l+ u k+l,l— v k,l+l~ v k+l,l) ; (3.7) 

(%-|-l,/-l-l 3^+1,/ yk,l-\-l ~Pyk,l) I 3^^,/— 1—1 3^+1,/^ 

= -^(Myt,/+1+ M fe+1,/+ V fe,/+1+ V fc+1,/) 

X («*+1,/+!H-Mjt+1 ,l——«*,/H-Vjfc+1,/+!H-Vit+1,/—Vjt,/+1 — V*-,/), (3.8) 


where 


(Zk+l,l+l -Zk+Ll ~Zk,l+l +Zk,l ) +Zk,l +1 — Zifc+1,/J 

= - — v*,/+t—vjt+i,/) 

x ( u k+l,l+l+ u k+l,l^ u k,l+l~ u k,l~ v k+l,l+l~ v k+l,l+ v k,l+l+ v k,l) • 


u k,l 


iln 


f fk,i8k,i \ 
\fk,igk,ij s ’ 


Vk,l 


iln 


f fk.lgk.l \ 
\fk,lgk,lj s ’ 


yk,i — ka— (In (fk,ifk,i))si z k ,i — ka — (In (g k ,igk,l))s, 


(3.9) 

(3.10) 

(3.11) 


Xk,i = yk,l+Zk,i = 2ka ~ (In {fk,lfk,lgk,lgk,l)) s • (3-12) 

The proof is presented in Appendix C. 

Finally we show that Eqs. (13.6b — (13.91) converge to the semi-discrete CSP equations (12.8b — 
(12.10b by taking a continuous limit in time (b —> 0). Under this limit, Eqs. (13.6b — (13.9b become 
d d 

-(uk+ 1 -Uk) + — (v £+1 -v k ) = (yjfc +1 -yk){u k +1 + ^ + ^+1 -hv*), (3.13) 

-j-(«jfc+l -Mjt) - ^r(v k+ i -V k ) = {Zk+I-Zk)(uk+1 + U k -V k+ 1 -V k ), (3.14) 

ds ds 
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^(yk+ 1 -yk) = ~{u k +i + u k + v k+ i +v k )(u k+ i-u k + v k+ i -Vjfc), (3.15) 
Y s ( Zk +1 = -^(Wfe+l+Mfe-V^+l -V fe )(Mfe + l-M^-Vfc + l +Vjt), (3.16) 


where Ff+ 2 fc Ff —y d s F(b —)• 0) is used. Obviously, we have 

-^-(Mjfc+i-Mjt) = -U*(m*+i + m*)-^( v *+i +v k )(z, k+ ] -z k -y k+] +y k )- J (3.17) 
ds 2 2 

^(v k+ i-v k ) = ^8 k (v k+ i+v k )-^(u k+ i+u k )(z k+ i-Zk-yk+i+y k ), (3.18) 


(3.19) 


from (13.131) — (13.141) . and 

~ x k ) = —^(“fc+t — + v fc+i “ v fc)’ 

by adding (13.151) and (13.161) . Eq. (13.191) coincides with Eq. (12.101) 

Finally, in view of the relations (IC.17I) - (IC.22I) . we have 

(z k +1 —z k —y k + i +yk)&k 
= ( z k+ i-z k - y k+ 1 +y k ) (z k + \-z k +y k+ \-y k ) 

' 2 r 'fk+lfk fk+lfk \ 2 f gk+\ 8 k gfc+lgAO 2 ' 

\fk+\ fk fk+lfk ' ygk+lgk gk+\gk' - 

= (uk+ 1 - U k ){v k+ 1 - v*). (3.20) 

A substitution of (13.201) into (13.171) — (13. 1 81) yields (12.8b — (12.91) . 

From the construction of the fully discrete analogue of the CSP equation, the multi- 
soliton solution can be expressed in the following determinant form 


a- 

T 


U k ,i = iln 

vjfc / = iln 


fk,igk,i 

fk,lgk,l 

fk,lgk,l 

fk,lgk,l 


i ( f k,I g k,l f k,I 8 k,l 


fk,l 8k,l fk,l 

,/ 


gk,l 


j_ / fk,l + 8 k,l _ f k,l _ 8'k,l 
2 l fk,l 8k,l fk,l 8k,l 


y k ,l — ha — (In (f k ,ifk,i)) s = ka — - 


Zk,i — ka — (In (g k ,ig k ,i))s = ka — - 



thus 


Here 


x k ,i = 2 ka - (In (fk,ifk,l8k,l8k,l)) s = 2ka-- 


1 ( f k,l 8 k,l f k,l 8 k,l 


■ + 


+ ■ 


+ ■ 


2 y fk,l 8k,1 fk,l 8k,l 

fk,l =^o(k,l) , fk,i = ii(k,l), gk,i = ^o(k,l), 8k,l = ^i(k,l), 

fk,i — Po(M ), /'*,/ = Pi (M), g' ktl = p' 0 (k,l) t g'k,i = Pi(k,l), 


(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 
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with x n (kj) and x f n (k,l) defined the same as (13.11) . p n (k,l) and p f n (k, l ) defined as 


P n(k,l) §( n+ j_ 2 ^(k,l) , P„(k,l) V(„+y-3) l) 


1 <i,j<N 

under the 2-reduction condition q t — —p\ (i = 1, ■ • • ,1V). 


1 <i,j<N 


(3.28) 


Remark 3.1. Two intermediate variables )y and Zk are introduced in constructing the fully 
discrete CSP equation. This often happens when we construct the fully discretizations of the 
coupled system such as the coupled modified KdV equation lf28l . 


4. Integrable self-adaptive moving mesh method 

In this section, we propose a self-adaptive moving mesh method for the CSP equation 
(1 1.2b — (11.3b and demonstrate the advantage of this integrable scheme by performing several 
numerical experiments. 

4.1. Numerical scheme 

One of the self-adaptive moving mesh methods for the coupled short pulse equation can be 
constructed by applying a semi-implicit Euler scheme to its integrable semi-discrete CSP 
equations (12.81) - (12.101) . The resulting numerical scheme reads 

p"C =pl+ y8?« + i + 4 )- + 4 ) , (4.1) 

«T‘ = 4 + ySjK-H + 4 ) - W+I + 4 ) . (4.2) 

5 f‘ = -J ((O 2 + (4ii) 2 - K +1 ) 2 - « +l > 2 ). (4.3) 

Here pk — Uk+i — Uk, qk = — vy, 8k — ay+i — ay- The superscript n represents the numerical 

value at t — nAl. The periodic boundary condition is applied. For convenience, we reserve 
the time t —* — t so that the left-moving wave becomes right-moving one. In what follows, we 
report the numerical results for one- and two- soli ton solutions. 

4.2. Numerical experiments 

For the sake of numerical experiments, we list exact one- and two- soliton solutions for the 
continuous, semi- and full-discrete CSP equation. 

(1). One soliton solution: the x-functions for one soliton solution of the CSP equation 
(11.21) — (11.31) are 

/ocl+ie 6 gocl+Fie 9 , 


(4.4) 
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where si = (1 — pi)/(l +Pi), 0 = Piy + s/pi + yo- This leads to one-soliton solution in 
parametric form 


u{y,s) — — (sech0 + sgn(si)sech(0-A)), 

Pi 

(4.5) 

v(y,s) — — (sechO —sgn(si)sech(0 —A)), 

Pi 

(4.6) 

x — y —-(tanh(0) + tanh(0-A)) . 

(4.7) 


Pi 


where exp(—A) = |si|. 

For the semi-discrete CSP equation, the x-functions are 

1 +ap l ^ p s/p,+y 0 


/toe l+i -- gS/Pi+yo gk od+i 5l 

-apij 


1 + ap | 


1 — ap | 

Finally for the fully discrete CSP equation, the x-functions are 


„s/p i+yo 


fu l+i 


1 +api\* /1 +hpj 1 


1 —ap\ J yi-bp^ 


(4.8) 


(4.9) 


8k,i x 1+isi 


/ l+api 
\l-api 


l+fr/>i 1 

1 -W 


l 


(4.10) 


The initial conditions for one-soliton propagation are taken from (14.41) with parameters 
yo — 0 , and p\ = 0.9, p\ = 2.0. The initial profiles are shown in Fig. [Ha) and (b), respectively. 
The simulations are run on a domain [—40,40] with 800 grid points, thus the average mesh 
size is 0.1. 

When p\ — 0.9, u is symmetric with two-spikes, and v is antisymmetric. The comparison 
between the numerical and analytical results is shown in Fig. [2] together with the nonuniform 
mesh. It can be seen that the non-uniform mesh is dense around the peak points of solitons. 
Moreover, the most dense part of the non-uniform mesh moves along with the peak point. 
When pi = 2.0, u is antisymmetric, and v is symmetric with a loop structure. The comparison 
between the numerical and analytical results is shown in Fig. [3] The error between the 
numerical solution and the analytical one is displayed in Fig. [4] 

(2). Two soliton solution: the x-functions for two soliton solution of the CSP equation 

(OHIO) are 


foe 1 +ie 01 +ie 02 — Z? 12 e 01+G2 , 
g 1 -\-is\e Qx -t-is^e 02 — h^e 01 ' 1-02 , 


with Si = (l-pi) / (l + pi), 0/ = piy + s/p i + }’; 0 (i = 1,2), andh 12 = (pi 
and b' n = bu*s\ *S 2 - 

For the semi-discrete CSP equation, the x-functions are 

k /1 i \ k 


fk 


l +i fi±f£!) e^ +i l±m 


1 — api 


\-ap 2 


Xi 


(4.11) 

(4.12) 

P 2 ) 2 /{pi + P 2) 2 ■> 












u(x,t) u(x,0), v(x,0) 
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X 


Figure 1 . Initial conditions for CSP equation, (a) p\ = 0.9; (b) p\ = 2.0. 




Figure 2. Comparison between numerical and analytical solutions for one-soliton to the CSP 
equation with p\ = 0.9 at t = 4.0; solid line: analytical solution, blue dot: numerical solution, 
red dot: self-adaptive mesh;(a) profile of u, (b) profile of v. 


8k' 


_ ( i}+ a Pl)^+ a P2) \ k b Mi 

l+isi (l±w±) k +iS2 


1 —ap i ) 




1 —ap2 J 


-b'v 


(1 + ap i) (1 + apf) \ k ? 1+ ^ 2 


12 ' (1 -ap\){\-apf)) 
with ^i = s/pi+yt o (i= 1,2). 

For the fully discrete CSP equation, the x-functions are 


fk,i « 1+ i 


1 +api\ k f 1 +bp x 1 


l-api) ll -b Pl l 


+ i 


l+ap2\ fl+bp 2 


1 -ap 2 ) \1 -bp 7 


(4.13) 


(4.14) 















































Error for u u(x,t) 
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Figure 3. Comparison between numerical and analytical results for one-soliton solution to the 
CSP equation with p\ = 2.0 at t = 12.0; solid line: analytical solution, blue dot: numerical 
solution, red dot: self-adaptive mesh; (a) profile of u, (b) profile of v. 
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Figure 4. The error between numerical and analytical results for one-soliton solution to the 
CSP equation with pi = 2.0 at t = 12.0; (a) error in u, (b) error in v. 


gk,! 0 * 


— b 12 

1+isi 


( (1 +ap\)(\ +api) \ k ( (l+bp l l )(l+bp 2 ! ) \ 
\(l-api){l-ap 2 )J y(\-bp- l )(l-bpf 1 )) 

/ I +apA k f 1 \ / 1 +ap 2 \ k ( 1 +bpf l 

\l-apij yi -bpi 1 ) lS2 \l-ap 2 ) yi -bpf 1 

( (l+api)(l+ap2) \ k ( {l+bpf l ){l+bp2 - 1 ) \ 
\{l-api)(l-ap 2 )J y(\-bp- l )(\-bpf 1 )) 


(4.15) 


l 


(4.16) 


As pointed in [fl6l . when p\ and p 2 are complex conjugate to each other, two-soliton 
solution becomes a breather solution. Eqs. (14.1 1I) - (I4. 121) are used as initial condition with 
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parameters chosen as p j = 0.4+ i, P 2 — 0.4 — i, y 10 = y 2 o = 0. The numerical results at t — 10 
are displayed in Fig. [5] in compared with analytical solution. Here a grid points of 800 is 
used on a domain [—40,40], the time step size is taken as At = 0.005. The error between 
the numerical solution and the analytical one is displayed in Fig. [6] It can be seen that the 
numerical results are in good agreement with analytical ones. 




Figure 5. Comparison between numerical and analytical results for breather solution to the 
CSP equation for p\ = 0.4+ i, p 2 = 0.4 —i at t = 10; solid line: analytical solution, dashed 
line: numerical solution; (a) profile of u, (b) profile of v. 
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Figure 6. The error between numerical and analytical results for one-breather solution to the 
CSP equation at t — 10.0; (a) error in m, (b) error in v. 
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In this paper, we proposed integrable semi-discrete and fully discrete analogues of a coupled 
short pulse equation. The determinant formulas of iV-soliton solutions for the semi-discrete 
and fully discrete analogues of the CSP equations are also presented. In the continuous limit, 
the fully discrete CSP equation converges to the semi-discrete CSP equation, then further 
converges to the continuous CSP equation. 

In a series of papers by one of the authors, we have constructed integrable discretizations 
for a class of soliton equations with hodograph transformation, and successfully used them 
as self-adaptive moving mesh methods for the Camassa-Holm equation [f30L [31 1 and the short 
pulse equation lfl6l . Based on the semi-discrete CSP equation (II.2b — (II.3b . a self-adaptive 
moving mesh method is constructed and used for the numerical simulations of the CSP 
equation. It should be pointed out that the feature of self-adaptivity of the mesh is due to the 
hodograph transformation. In other words, the hodograph transformation converts the uniform 
and time-independent mesh into a non-uniform and time-dependent mesh. It is a further topic 
to seek for such kind of self-adaptive moving method when the hodograph transformation is 
not present. The numerical results confirms that it is an excellent scheme due to its nature 
of integrability and self-adaptivity of the mesh. This is our first time to extend this superior 
numerical method to a coupled system. 
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Appendix A. Proof of Proposition 1 

Proof. For simplicity, we introduce a convenient notation 




§n\k) ^i+jV-lW 

$n\k) tfl(k) ••• ^_t(*) 

^\k) ••• 1 (*) 

vi 1} (£) vi+i(£) vi+iv-i(fc) 
,J 2 )tlA ™( 2 ) dA ... „,( 2 ) (lA 


(A.l) 


(A.2) 


Since 


d X - l ^(k) = <bnU’<)i <bn\k+l)-$\k) = a$ +l (k+l), 
dx_! (k) = vjlj (k), (k + 1) - (k) = avS-t ( k + 1), 


(A.3) 


(A.4) 
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${k)~vlP{k)=$® 1 (k), (A.5) 

we can verify the following relations 

d x _Mk) = I —U, h, ■ ■ ■ ,N-2 k ,N- h\ , (A.6) 

t n +i(Hl) = \h,2k,---,N-2 k ,N-lk,N k + 1| 

= 1 14,2jt, • • ■ ,N — 2 k ,N — l k ,N — l k+i | , (A.7) 

T n (k+1) = |0jt+i, 1 at+i, - - - ,7V — 2^ + i,A^— ljfc+i| 

= |0jfc, Ik, - ■ ■ ,N — 2k,N — ljt+il , (A.8) 


d x _ 1 'Zn(k+ 1 ) = |—lfc, lyt, •• ■ ,N — 2 k ,N — lfe + i| + | 0 jk, l k , - ■ ■ ,N — 2 k ,N — 2 k+l \ 

= | —U, U, •• ■ ,N — 2 k ,N — lfc + i| +a\0k, U, •• -,N — 2k,N— ljt+i| ■ (A.9) 

Combining (IA.8I) with (IA.9I) . we have 

( - a*.! - A x„ (* + 1) = - I - h, u, • • •, N - 2,, N - 1* +1 1 . (A. 10) 

\a J a 

Therefore, the Pliicker relation for determinants 

|0*, l kl ---,N-2 k ,N-l k \\-l k ,l k ,---,N-2 k ,N- U+i | 

+ |U,---,A-2,,A-U,A-U +1 ||-U,0 fc ,U,---,A-2,|=0 (A.ll) 

gives 

(-d x _! ~ 1 ) i n (k+ 1) X x„(fc) - -%n(k+ 1) X a x _jX n (k) +x n+ i(k+ l)x„_i(k)-0, (A. 12) 
\a J a 

which is nothing but the bilinear equation (12.11) . Eq. (12.21) can be proved in the same way. 
Now we proceed to the proof of Eq. (12.31) . Similarly we can verify the following relations 


x n(k)= \0 k ,l k ,---,N-2 h ,N-l k \ 


1 

<N 

i 

O 

(A.13) 

x»+t(*) = \l' k ,2' k ,-,N-l' k ,N k \ 


1 

1 

(N 

(A. 14) 

(3^,-1 )x„W = |-li, . 

(A. 15) 
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Therefore the Pliicker relation for determinants: 

\-\' k , \' k , ■ ■ ■ ,N-2 k ,N — l k \ \0’ k , l k ,---,N — 2' k ,N - l' k \ 

- \0 k ,l k ,---,N-2' k ,N-l k \ \ — \ k ,l k ,---,N — 2 k ,N—l' k \ 
+ \-ll0lll...,N-2[\\ll---,N-l' k ,N-h\=0, 


gives 


(9jc_! - 1 )*«(*) X <(fc) -Tn{k) X d x _X(k) +T n+ l(k)f n _ l (k) = 0. 
Therefore, Eq. (12.3l) is proved. 


(A. 16) 

(A.17) 
H 


Appendix B. Proof of Theorem 1 

Proof. By putting s = 2x-i, %o(k) — f k , Xi (k) — f k , (12.11) can be converted into 


D s - 1^ fk +1 ■ fk = -fk+ifk, 

-D, s -\)f k+] -f k = -f k+ vfk, 

a J 


-D s — 1 gk +1' gk = -gk+igk, 

1 J 
\ 

D s — 1 g k+ 1 • g k = . 


(B.l) 

(B.2) 


while by putting %' 0 (k) — g k , ( k ) = g k , (12.21) can be converted into 

'2 
—j 
a 

'2 . 
a 

Furthermore, the above bilinear equations can be rewritten as the following logarithmic 
derivatives 

2 


(B.3) 

(B.4) 


ln fk+^ 

fk j 

l -- 

A+t/jt 
fk+lfk ’ 

(B.5) 

In 

. fk , 


A+t/it 

A+t/fe ’ 

(B.6) 

In 8k+] 

. gk y 

).- 1 = ’ 

gfc+tgfc 
g/fc+lg/t ’ 

(B.7) 

In 8k+] 

. gk y 



(B.8) 


Introducing two intermediate variable transformations 
<5 k {s) = 2iln , o' k {s) = 2iln 

one arrives at a pair of semi-discrete sine-Gordon equations 

1 , N / a* +) +aA 

— (<J»+1 -o k ), = sin I - j , 


( gk(s) 

\gk{s) 


(B.9) 
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It then follows that 



where c^ iS denoted the derivative of a/ t - with respective to s. 
Next, we introduce dependent variable transformations 



and discrete hodograph transformation 

x k = 2 ka - (In {fkfkgkgk)) s ■ 

Then the nonuniform mesh can be derived as 


= 2 a 


S/t — ^-k+ 1 

fk+ lfk+lgk+lgk+1 \ 
fkfkgkgk J ■ 


In- 


a_ f fk+lfk fk+lfk gA+lgA gfc+lgA 

2 V/a+i/a /a+i/a Ia+iIa ?4+l?t/ 

_/ 


= a cos 


/ <*A+1 + &k 

V 2 


+ cos 


<%1+ C A 


Taking the derivative with respect to s results in 

dh ' / - 1 - x 




= a 


^cos ^ 


GAM-! +Gk 


T cos 


°A+l+ a A 


= -^(4 +l -< + v 2 k+l -v 2 k ). 
Furthermore, assuming 

'ga + Oa+i + o'+g' + i 


p k = sec 


<?A = sec 


^A + ^A+l “G, 


J A+1 


(B.10) 

(B.ll) 
(B. 12) 

(B.13) 

(B.14) 

(B.15) 


(B. 16) 


(B. 17) 


we have 

™ 2 a 

Ok = -, 

Pkdk 


(B.18) 
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dp , 1 d (o k + a k+l +a' k + a' k 


— — cos 


ds ds 


k +1 


. f + °£+i + o' k + o[ j \ U k + U k+ 1 

— sin ---— --- 




dk_ 

2 


sin 


f Gfc + Gfc+i \ 






sin 


+ ®k +1 \ \ u k + u k+ \ 


dk 1 / | / / \ + 

- T ^ (a * +1 ~ a * + ° fc+1 " 


c lk / 2 


4a 


■ 4 ), 


and similarly 


dq k 1 


P* / 2 




)• 


ds 4a 
Thus, in turn, Eqs. (IB. 19b and (IB.201) become 


d PZ 2 u k+i- u l 


and 


ds 

dq\ 


= Pk- 


8/t 


(B. 19) 


(B.20) 


(B.21) 


= d k 


2 V k+l V k 


ds 1K 8 k 

respectively. 

On the other hand, with the help of trigonometric identity, p k can be expressed as 

2 f O k + O k+l + c' k + o[. + 1 


(B.22) 


p k = 1 + tan" 


= 1 + p% sin 2 


c k + c k+ i+a' k + o' k 


k +1 


= 1 + 


phi (,.i 


4a 2 


hk+l ~ U k) 


= l+(^± 


-U k 


\ h 


and similarly 

dk = 1 + 

Therefore, we finally have 

d (u k +i ti k 
ds \ 8jt 

d ( v k+l — V* 


Vjt+l ~Vk 


ds 




= 1 + 


= 1 + 


^ Uk +1 | Mfc+1 +U k 


V k +\ ~ v k 


V k +1 +Vjfc 


(B.23) 


(B.24) 


(B.25) 


(B.26) 


Substituting (IB. 171) into (IB.25I) - (IB.26I) . one arrives at (12.8b — (12.9b . the first two equations of 
the semi-discrete CSP equation. H 
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Proof. The bilinear equations (13.2b — (13.51) imply the following bilinear equations 
{z Ds ~ ^ fk+\,rfk,i + fk+\,ifk,i = 0, 

(~ Ds ^ ^ fk+1,1 ' fk,l + fk+1,1 fk,l = 0, 

gk+i,rgk,i + gk+u8k,i = 0, 

(~ Ds ~ ^ 8k+l,l'8k,l + 8k+l,l8k,l — 0, 

(2 bD s - l)fk,i+i ■ fk.l + .fk.lfk,l+ 1 =0, 

(2 bD s - l)fk,i+i ■fu + fk.lfu+\ — 0, 

(2 bD s - l)sk,i+i-8k,l+8k,l8k,l+i — 0? 

(2bD s - l)gpi +i ■ g k j + g kJ g k) / + i = 0, 

which can be rewritten by logarithmic derivatives as 


2 

a 

2 

a 

2 

a 

2 

a 

2b 

2b 

2b 


In 


fk+1,1 

fk,l 


In -fk+\-l 


= 1 - 


- In 


In 


fkj / a 
ft+u' 

8k,l 
8k+l,l 
8k,l 


= 1 


= 1 


( ln fuu 

V fk,l 

(i„C±i 

V fk,l 
( in^Mti 

8k,l 


fk+l,lfk,l 
fk+l,lfk,l ’ 
fk+l,lfk,l 
fk+l,lfk,l ’ 
gfc+1,/£*:,/ 
a 8k+\,l8k,l 

_ j _ gfc+Ug/M 

fk,lfk,l +1 
fk,l+\fk,l ’ 
_ _ fk,lfk,l+l 
s fk,t+lfk,l 


1 - 


26 (In 


= 1 




V 8k,i 


gfe,/gfc,/+l 

s 8k,l+l8k,l 

_ | _ £fc,/gfc,/+l 

a 8k,l+l8k,l 


(C.i) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 

(C.7) 

(C.8) 

(C.9) 

(C.10) 

(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 

(C.16) 


Based on the dependent variable transformation (13.101) and discrete hodograph 
transformation (13.1 II) . we can verify the following relations 

fr , f, , a, ., , a, ,N 

(C.17) 
(C.I 8) 
(C.19) 


Uk+l,l~ u k,l — 

i a 

( fk+1,1 fk,l 

, 8k+l,l8k,l 

fk+1,1 fk,l 

8k+l,l8k,l\ 

2 1 

{ fk+\,lfk,l 

1 

§k-\- 1 ,1 §k, l 

fk- (-1, / fk , / 

gk+\,lgk,l J 

Uk,l+\+Uk,l — 

i 

( fk,ifk,l+l 

, 8k, 18 k, 1+1 

/it,/ fk,l+ 1 

8k,l8k,l+ 1' 

2b 

^fk,l+lfk,l 

1 

8k,l+\8k,l 

fk,l+lfk,l 


Vk+l,l~Vk,l — 

i a / 

'fk+1,1 fk,l 

8k+l,l8k,l 

fk+1,1 fk,l 

j ^ 

2 \ 

' fk+ 1 ,1 fk, l 

8k+l,l8k,l 

fk+1,1 fk,l 

“ ) 

g/t+Ugfc,/' 
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i ,fk,lfk,l +1 8k,l8k,l +1 fk,lfk,l +1 , 

Vt /_i_i + Vt / = - -=-=-1- 

26 fk,l+lfk,l 8k, l+l 8k, l fk,l+lfk,l 8k,l+l 8k,l 

a (fk+l,l fk,l . fk+l,lfk,l 

yk+l,l~yk,l = 7T 1 --7-1- — 


2 'fk+l,lfk,l fk+l,lfk,l' 

a (8k+l,l8k,l , 8k+l,l8k,l\ 

Zk+\,1-Zk,l = ■= -h--—I, 

2 \gk+\,i8k,i 8k+i,i8k,r 

_ 1 1 ( 

yk,t+i-yk,i -~b + 2 b\ 


Zk,l +1 Z k ,l + 


1 , 1 ( fk,lfk,l +1 ^ fk,lfk,l +1 
'fk,l+lfk,l fk,l+lfk,l' 
1 . 1 (8k,l8k,l+\ 8k,18k,l+l 


+ ■ 


b 2 b \8k,l+l8k,l 8k,l+l 8k,l ‘ 

Then, the ratios on the r.h.s of Eqs. (1C.17MC.24I) can be solved as 


fk+l.lfk,l _ 1_ 
fk+l.lfk,l Cl l 
fk+l.lfk,l 1 


fk+ 1, / fk. I a 
fk , / yir, 1 


yk+ 1,1 -yk,l - 2 ( u k+l,l - u k,l+Vk+l,l —Vk,l) 

i , ' 

yk+ 1,1 —yk,l + +.\ u k+l,l - Uk,l+ v k+l,l ~ v k,l) 


fk,l+ \ fk,l 
.fk, l fk,l +1 


= 1+6 

= 1+6 


fk,l+lfk,l 

8k+l,l8k,l _ -j_ 
8k+l,l8k,l Cl L 

_ -l 

8k+i,i8k,i a L 
8k,18k,l+l 


yk,l+ 1 - )+/ - 2 ( u k,l+l + u k,l + Vk,l+ 1 + Vfc,/) 
3+/+1 — Jyt,/ + 7 :( m A',/+1 + «£,/ +Vk,l+1 +Vk,l ) 


Zk+l,l~Zk,l - 7t( u k+l,l - Uk,l —Vk+l,l+Vk : l) 


Zk+l,l~Zk,l + 2^ Uk+1 ’ 1 ~ u k,l ~ v k+l.l +Vk,l) 


8k, l+l 8k, l 
8 k, 18 k, l+l 


= 1+6 
= 1+6 


Zjfc,/+1 _ z k,l ~ o ( u k,l+ 1 + Uk,l ~ Vk,l+ 1 - Vjt,/) 


Zjfc,Z+l - Z/t,/ + ^ (+M+1 + Uk,l ~ v k,l+1 ~ V k ,l) 
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(C.20) 

(C.21) 

(C.22) 

(C.23) 

(C.24) 

(C.25) 

(C.26) 

(C.27) 

(C.28) 

(C.29) 

(C.30) 

(C.31) 

(C.32) 


8k,l+l 8k,l 

By making a shift of /—>■/+ 1 in (IC.25I) . then dividing it by (IC.26I) . meanwhile, dividing 
(IC.27I) by (IC.28I) after a shift of k — > k + 1, one obtains 

3+H./+1 —yjt,/+l - 4(+t+l,/+l — “£,/+! + V/fc+l,J+l “ v E/+l) 


1+6 


Jyt+l,/ - yk,l + ^{tlk+ 1,1 - U kt i + Vk+ 1,1 - v k j ) 

yk+l,l+l —yk+ 1,1 — j{Uk+l,l+l +Uk+l,l + Vk+l,l+l +Vk+l,l) 


l+b yk,i +1 — yfc,/ + 5 (Mjt,/+i + M fc,z +v&,z +1 + +t,/) 

Similarly, one can obtain 

z*+i,f+i ~Zfe,/+i - i(++-U+i v^+i,/+i + vfc,/+i) 

Zk+ 1,1 - Zk,l + j{u k +l.l - u k j - V k+ 1 J+V k j) 


(C.33) 


1+6 

Zk+l,l+l —Z k + 1 ,/ — j(u k +l.l+l+U k +ij —V k +iJ+i — 

n+i,i) 

1+6 

Zk,l+ 1 -Zfe,/ + 5 («*,/+! + +t,/ - V ^/ + 1 - Vjt,z) 



(C.34) 
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from relations (1C.29I) - (1C.31I) . Equating the real parts and imaginary parts of (IC.33I) and (IC.34I) . 
we have the fully discrete CSP equations (13.61) — (13.9b . H 
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